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Abstract 

At sufficiently high baryon densities, the physics of a dense quark- gluon 
plasma may be investigated through the tools of perturbative QCD. This 
approach has recently been successfully applied to the study of color super- 
conductivity, where the dominant di-quark pairing interaction arises from one 
gluon exchange. Screening in the plasma leads to novel behaviour, including 
a remarkable non-BCS scaling of Tq, the transition temperature to the color 
superconducting phase. Radiative corrections to one gluon exchange were pre- 
viously considered and found to affect Tq. In particular, the quark self-energy 
in a plasma leads to non- Fermi liquid behaviour and suppresses Tq- However, 
at the same time, the quark-gluon vertex was shown not to modify the result 
at leading order. This dichotomy between the effects of the radiative correc- 
tions at first appears rather surprising, as the BRST identity connects the 
self-energy to the vertex corrections. Nevertheless, as we demonstrate, there 
is in fact no contradiction with the BRST identity, at least to leading log or- 
der. This clarifies some of the previous statements on the importance of the 
higher order corrections to the determination of Tq and the zero temperature 
gap in color superconductivity. 

PACS numbers: 12.38Aw, 12.38.-t, ll.10.Wx, 11.15.Ex 
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I. INTRODUCTION. 



Attention to the physics of a dense quark-gluon plasma has recently been revived, along 
with progress in understanding the phase structure of QCD. Interest has only heightened 
recently with the projected onset of relativistic heavy ion collisions at BNL. Such unusual 
conditions of QCD may also exist in nature, for example in the core of a dense neutron star. 
The asymptotic freedom of QCD makes a perturbative treatment applicable at sufficiently 
high baryon density, and the attractive di-quark interaction mediated by one gluon ex- 
change in the antisymmetric color representation induces superconductivity below a certain 
temperature [|TJ— 0] . 

Working at non-zero temperature and chemical potential introduces several complica- 
tions. One of the primary features of the plasma is that it screens the QCD interaction. 
Thus it is necessary in principle to dress gluon propagators with hard dense/thermal loops 
(HDL/HTL) in the plasma. For conditions in the range of interest for color superconduc- 
tivity, the temperature effects are less important, and only the effects of screening by HDL 
need to be considered. That this screening is important is a corollary of a more general 
statement, namely that, with a long range interaction at non-zero chemical potential, a 
straightforward power series expansion of the free energy in the coupling g results in infra- 
red divergences. A resummation over the fermion loops, and the replacement of the bare 
gluon propagator by that dressed with HDL [fl-|T0|l, prior to perturbative expansion resolves 
the infra-red difficulties. This was demonstrated, for example, for a non-relativistic electron 



gas with Coulomb interaction in [11]. The resultant perturbative series contains logarithms 
of the coupling constant accompanying the powers of it. 

As a result of HDL, the electric gluon propagator is screened effectively by a Debye mass, 
mo-, while the magnetic propagator is poorly screened via Landau damping in the sensitive 
region of momentum space. An important consequence of this is the introduction of non- 



Fermi liquid behavior of the quark self-energy [12]. In the infra-red limit (highlighted by a 
cutoff l c ) this self-energy is 

ia 2 4/ 3 
S (<^U--I^7 4 ,ln^, (1.1) 

where (p, v) are the Euclidean energy-momentum. Such a non-analytic dependence on energy 
v was first discovered in solid state physics [13] in the context of magnetic interactions. The 



logarithm suppresses the quasi-particle weight at the Fermi level and the single fermion 
occupation number becomes a continuous function at p = pp, the Fermi momentum, in 
contrast to the kink of a Fermi liquid. Another effect is a term ~ TlnT in the specific heat, 
but this turns out to be too small to be observed since the magnetic coupling represents 
merely a relativistic correction. Such non-analyticity as indicated by (|1.1|) was also suggested 
for certain strongly correlated systems such as high Tq superconductors [pj . 



In a relativistic quark-gluon plasma, the relatively poor screening by Landau damping 
is far more transparent. The di-quark pairing force is dominated by magnetic gluons, and 
Landau damping gives rise to a remarkable non-BCS scaling of the transition temperature 
and the energy gap for color superconductivity |l2| , |l5| -p^ , 



k B T c = c^e V ^ o . (1.2) 

sr 
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The non-Fermi liquid behavior (|1.1|) suppresses the pre-exponential factor c significantly, 
and we found that [ lq 



Coe -^(- 2+4 )(^- 1 ) ~ 0.176 c for N c = 3, (1.3) 



with Co the pre-exponential factor without radiative corrections |lq , P^Jll^ , |20[| . We have also 
argued that the contributions from other higher order diagrams to ( |1.3[ ) are subleading. 

Since the inverse quark propagator is related to the quark-gluon vertex function through 
a BRST identity, so is the quark self-energy £(P) of Fig. |I| to the radiative correction 
AJP',P) of Fig. |2|. To illustrate such a relation, we focus on the vertex correction in 
Fig. ^|a, which survives in the abelian case and will be referred to as the 'abelian' vertex 
in the following. After factorizing out the group theoretic coefficients from the vertex and 
self-energy, 

Af\P', P) = g TfT l f TfA^P', P) (1.4) 

and 

E(P)=T l f T l f E(P), (1.5) 

we have the Takahashi identity 

(P'-P),A,(P',P) = E(P')-E(P). (1.6) 
Taking the limit P' — > P, we end up with the usual Ward identity, 

A /I (P,P) = ^-H(P), (1.7) 

which is similar to the Ward identity of QED. Because of the above behavior of S(P), this 
identity raises a suspicion that A At (P', P) must also contribute to the pre-exponential factor 
in a manner similar to ( |1.3| ). This was, however, ruled out in [IS] where we showed that 



while the derivative J^S(P) contains the logarithm of u, the derivative Tp-S(P) does not, and 
thus the effect is not that of wave function renormalization. The contribution of A 4 (P',P) 
remains subleading even with the logarithm since the Coulomb propagator attached to it is 
strongly screened. 

Nevertheless a paradox arises here. The integral representations of A;(P, P) and A 4 (P, P) 



look identical while the above results, together with (|1.7|) , suggests different answers. In this 
article we shall disentangle this mystery. It turns out that the expression A M (P, P) is highly 
ambiguous in the presence of a Fermi sea, and in particular, 

lim lim AJP',P) + lim lim A^(P',P), (1.8) 

p'—*pv'—*v u'^up'^p 

a common ambiguity of the infra-red limit (the zero energy-momentum limit of soft lines 
of a diagram) in the absence of covariance. The contribution of A^P', P) to color super- 
conductivity comes mainly in the region \v' — v\ <C \i with \p' — p \ ~ m^ 3 |z/ — z^ 1 / 3 while 
\p — n\ ~ \v\ and \p' — ji\ ~ \u'\, which is closer to the order of the left hand side of ( p..8| ). 
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FIG. 1. The quark self-energy diagram. 

By carefully tracing the subtleties of the infra-red limit along the different routes, we are 
able to reconcile the logarithmic behavior of (|1 . 1| ) with the Ward identity (|1.7| ) as well as 
the full BRST identity when the group theoretic factors and the vertex diagrams in Fig. 
and Fig. |]c are restored. Yet the suppression in (|1.3|) remains intact. 

Though we are mainly addressing QCD and color superconductivity in this article, the 
non-Fermi liquid behavior of the fermion self-energy and the vertex function apply, to a 
simpler extent, to the relativistic electron plasma as well. Such a plasma exists inside a 
white dwarf star, a supernova or a red giant star, for which the condition that the chemical 
potential is much higher that the temperature is valid. 

In the next section, we shall calculate the quark self-energy and pin down the mathemat- 
ical mechanism behind the logarithm of The vertex function A 1 (P 1 , P) is analyzed in 
section III in light of the BRST identity. The contribution to color superconductivity will 
be discussed in section IV together with some concluding remarks. 



II. THE QUARK SELF-ENERGY. 

Little motivation is required for the analysis of the quark self-energy, represented in 
Fig. |l|; it may appear as a simple radiative correction in perturbative processes, but it also 
enters the BRST identity. The form of the self-energy also characterises the non-Fermi 
liquid behaviour at high density and has a subtle influence over the divergences of the quark 
vertices. Without more ado, in Euclidean space we write, 

1 r d 3 T 

S(P) = -g 2 T)T)- 53 J j^y 3 ^(Lh,S(L + P) 7 „, (2.1) 



where L = (/, —u n ), P = (p, —v n ), w„ = ne, v n — ( n + |J e, e = 2nk B T and TjTj = Cf 

"^T 1 for Tf in the fundamental representation of SU(N C ). Following the notation of [18,20 
we write the quark propagator as 

S(P) = L (2-2) 

where W = 74 (/i + iv) — ij • p. In the presence of a Fermi sea it is necessary to incorporate 
HDL into the gluon propagator at leading order. While it is possible that a magnetic mass 
of order T exists, at high density fi ^> ksT the damping due to HDL prevails over that due 



to HTL |2T|. Incorporating HDL, the gluon propagator in the covariant gauge takes the 
form 

V ^ K > = K 2 + a M (k,u) P ^ + K 2 [l + a E {k,uj)jk T \ >ilv ~ ^JlOf' (2 ' 3) 
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where K = (k, -u), K 2 = k 2 + u 2 , P? = 5 l3 - k t k v P£ = P£ = Pj A = 0. 



, ^ Kjjty _ pT 



(2.4) 



and a is the gauge parameter (we have adopted the notation \k\ = k). The electric self- 
energy a E (k,u) and the magnetic self-energy a M (k,u>) in (|2.3| ) are given by a E (k,u) = 
m 2 D f E (u/k), u M (k, u)) = m 2 D f M (w/jfe), with m 2 D ~ and 



/ M (*) 



1 — x tan 1 



x 



x 



1 + x 2 ) tan 1 



— x 



(2.5) 
(2.6) 



For more discussion of our notation or HDL in general see |T8| , pO| or [plj] , respectively. From 
( j2.5|) we can see that the Coulomb interaction is strongly screened while from (|2.6| ) the 
magnetic interaction is not. In this paper, we shall be interested only in the leading infra- 
red behavior, which comes solely from magnetic gluon exchange, and so we shall neglect the 
electric contributions and regard 



with 



k 2 + u 2 + a M (k,u)' 



(2.7) 



(2.8) 



To focus upon the infra-red behavior, we separate the loop integral into two regions and 
rewrite the self-energy as 



Z(P) = C f [E<(P)+E>(P)} 



(2.9) 



where the superscripts denote integration inside and outside the infra-red sensitive region: 
< I < l c , —ijj c <ijj n <uj c with l c ,u c <C [i. We shall evaluate the infra-red sensitive region 
only, 



E<(P) 



An 2 Jo 



lc 



dll 2 d(cos6)-J2 



74 - i(l -P) 7 -V 



(2.10) 



where £ = \l + p | — //. Corrections due to the change from a discrete sum to an integral 
are sub-leading; they can be obtained using zeta-function techniques as demonstrated for 
similar processes in [[HJ . Thus, in pursuit of the leading order behavior only, we immediately 
move to the continuous energy limit with k^T <C v m <C uo c . Making the change of variables 
from 9 to £, 



d(cos 0)[74 — i(l ■ p) 2 7 • p] 



lA-i-^a-ipfl -p 
l z p z 



(2.11) 
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with £ p = p — /J,, we find, 
E < (u,p) 



F(u,p;l,u) 



8tt 3 



dll / dujV(l, u)F(u,p; Z, to), 



d ^ FY , N • 



(2.12) 
(2.13) 



Fixing p = [i for the external lines and carrying out the integration over £, we have, 

Z \ 



7-,/ 7 \ CS- -1 ^ ^ C<J + Z/ ^ - / -, + Z/ _] 

i* (z/, //; I, uj) = 2^74 tan — ; h 2 — ; — 7 • p ( 1 : — tan 



UJ + V 



I 



I 



U) + V 



(2-14) 



so that 



— F(z/, /i; Z, = 2^7 4 (5(z/ + w) - _^ i/>2 ^ /2 : 1 



+ 



(w + v) 2 + Z 2 

^ + z^ _j I 2{u + v) 2 + l 2 
—2 — ; — tan h 



l 



1-p, 



(2.15) 



where the delta function comes from the discontinuity of the inverse tangent function. We 
find the energy dependence of the self-energy by differentiating, 



d_ 

dv 



E < (iy,p) 



= g 2 [A(v) + B(v)}, 



p=/i 



with 



My) = -7-574/ dl 



4tt 2 Jo p + u 2 + m 2 D fM (=vj 



(2.16) 



(2-17) 



and 



Bly) 



47T 3 



dl / dwV(l,uj) 



(2.18) 



il 



(u + v) 2 + / 



2 74 



'a; + 1/ 1 / 2{u + v) 2 + l 2 
2 — : — tan" 1 h 



l 



uj + u (co + u) 2 + l 2 



Noting the asymptotic behavior, that / (x) ~ | |x| for |x| <C 1, a scale /o may be introduced 

to divide the integration in A(v) into two: \v\ <C Z <C (m^z/l) 1 ^ 3 . For Z < Z we have the 
contribution, 



Wo 



Z 



< 



lo 



dl- 



l 



< 



1 



lo 



dll 



/ 3 



m D \v\ 



< 1. 



(2.19) 



All of the inequalities follow straightforwardly from the definition of Iq except for the second, 
which is due to f M {y /Z) being a monotonically decreasing function of I. Therefore, neglecting 
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this subleading contribution, we find a logarithmic infra-red singularity in A(u) arising from 
the second region (namely the region l < I < l c ) . The integration B{y) is finite in the limit 
v — > 0. We end up with 



4tt 2 



74 



dl 



I 2 



lo 



I 3 + vH + $m 



D\ 



12n 



r74 In • 



4/? 



+ ■ 



[2.20) 



That the self-energy does not depend upon the spatial momentum in the infra-red limit 
can easily be ascertained by differentiating H(P) with respect to pi. Noting that 
from ( |2.12|) and ( [2.13| ) it is straightforward to find, 



d_ 



ig 2 B(v), 



(2.21) 



p=/j, 



which is both real and finite in the limit v — > 0. Therefore the logarithmic singularity can 
not be attributed to a wavefunction renormalization. This is also the case found in a solid 



state physics context 13 



To summarise, we find that in a dense quark-gluon plasma the quark self-energy exhibits 
non- analytic behavior only for the energy component, 



E(P)| 



p=ii 



12vr 2 



C/74 i^ln 



4/ c 3 



+ ■ ■ 



(2.22) 



(the imaginary part of the self-energy, contributing to damping in the plasma [p2p3| , is 
analytic as v — ► 0). It is important to note that the infra-red non-analyticity in the energy 
originates in the discontinuity of the pole cutting the contour in the £ integration. This 
feature gives rise to the ^-function in (|2.15|) which ultimately leads to the infra-red non- 
analyticity. This behavior will be seen to repeat itself in section IIIA where it will lead to 
infra-red divergences in the radiative corrections to the quark-gluon vertex. 

From the result Q2.22p it is clear that covariance is broken; this is a direct effect of the 
presence of a Fermi sea. We also see that S<(0) = 0, so that the self-energy leads to no 
chemical potential renormalisation from the infra-red side. What is also of considerable 
interest is the BRST identity. How this is met in the dense quark-gluon plasma is subtle 
and we investigate this phenomenon in the next section. 



III. THE BRST IDENTITY AT HIGH DENSITY. 

Since the quark self-energy is only non-analytic in the external energy, one may expect 
from a generalisation of the Ward identity of QED that the Coulomb-quark vertex has 
similar behavior and is divergent while the magnetic gluon-quark vertex is not. However, 
it is also clear that the integrands in the infra-red region for and A- are mathematically 
identical, save for an indiced prefactor. Apparently we have something of a paradox: from 
the self-energy we expect only the Coulomb vertex to be divergent, but there appears to be 
no mathematical difference between the infra-red contribution to the Coulomb and magnetic 
vertices. We shall resolve this paradox in this section and show how the BRST identity works 
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at high density. First of all, we shall examine the zero energy-momentum transfer limit of 
the abelian vertex A jU (P', P) and derive the precise expression of the Ward identity (|Lj 



d 

lim lim A 4 (P',P) = -^r-S(P), (3.1) 

v'—>u p'^p ov 

d 

lim lim KdP'P) = -^-S(P). 



The latter relation for the magnetic vertex was previously investigated in |24] for a system 



of fermions interacting with transverse abelian gauge bosons. Contrasting the Coulomb 



and magnetic cases in (3J.) provides an important clue hinting that the ordering of limits 
contributes to the subtlety we are overlooking. To show how the paradox is resolved this 
identity shall be considered in the next subsection, where the infra-red behavior of the 
vertices is discussed and the abelian vertex is treated in detail. In the second subsection we 
shall show how the full BRST identity works in the dense quark-gluon plasma in terms of 
Feynman diagrams. 



A. Infra-red Behavior of the Abelian Vertex. 

To explore the behavior of the quark-gluon vertices and their relation to the BRST 
identity we shall analyse in detail the abelian vertex A^°) shown in Fig. 0a. We refer to 
this vertex as 'abelian' since it is the only physical vertex that also appears in the abelian 
theory. In order to simplify matters further, in this subsection we shall put both external 
quarks on-shell, p = p' = fi. 

Intuitively, we may expect the behavior of the vertex to depend subtly upon the ordering 
of the limits. We may develop this intuition from HDL, for example, where although there 
is an analytic result for the screening ( |2.5| ) and Q2.6Q , it has different asymptotic behavior in 
the two orderings of the limits (x —>■ and x — > oo). As we shall see, HDL and the BRST 
identity at high density are intimately connected and it is no surprise that the ordering of 
limits in (|3.1|) is crucial in resolving the paradox. We write the abelian vertex as 



Af\P',P) = gTpT l f TpA^{P',P) 

= sT;(-%^ + C7W)(P',P), (3.2) 



where 

A { ;\P',P) = l -jHj ^V vp {l,u) lv S{L + P') llx S{L + P) lp . (3.3) 

This may be written in terms of two integrals, one inside and one outside the infra-red 
sensitive region; < I < l c , —uj c < oj < u c with l c , oj c <^ /i, 

AW (P', P) = P M [A (a)< (P', P) + A (a)> (P', P)] . (3.4) 

with P^ = (— ip, 1). In this subsection we are only interested in the leading infra-red behav- 
ior. So we evaluate, 
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A (a)< (P',P) 
A(P,P';L) 



in 3 Jo 



dll 2 J d(cos0)(74 -t^ ■pcos 2 6)A(P,P';L), 



1 



duV(l,u) ( 1 



1 1 

+ 



Ib(-w) 



(3.5) 
(3.6) 



where £ = i/ + z£, £ = |Z + p \ — // and (' and £' refer to u', p'. The logarithm in ( p.6|) 
introduces a branch cut which we may take to lie along the positive real axis and the 
contour to run above and below it in the normal fashion. As shown for the self-energy, it is 
only the discontinuities that occur as poles cut the contour and branch cut that induce the 
infra-red singularity. Hence we need only focus upon the second and fourth terms in ( p.6|) , 
since the other terms are regular. Using the convention that arg(— 1) = 0, we find that the 
contribution of these poles reads 



A disc .(P',P; L) 



71 



c-c 



[Sign(0^(/,C)-Sign(O^(/,C')] 



(3.7) 



where the sign function comes from the discontinuity of ln(— u) crossing the cut. 

We are now in a position to take the limit P M — > P'^ for the Ward identity and we shall 
consider the two different orderings of the limits in turn: 







lim liniA disc .(P',P;L) = -7rSign(£)— C) 



u'—>u p'^p 



d 



lim. lim A disc .(P', P; L) = m— [Sign(0Z>(Z, C)] • 

" — " "■ — " at, 



(3.8) 
(3.9) 



In both cases we are looking at the infra-red limit, and thus fix the external momentum to 
be p = p' = fi, £ p = 0. 

First of all, considering case i), using the change of variables ( 2.11| ) it is straightforward 
to find, 



lim lim A (a)< (P',P) 

v'—*v p' '—>f> 



9" 



P=IM 



In 2 Jo 







dll I ^ S ign(0^P(/,C)(74-^7-p) 



W 
4tt 2 

ig 2 

YItx 



74 



r74 hi • 



dll [D(l, u) - V(l, v + + 
All 

-5 e r-r + "- J 



7xm D \i> 



(3.10) 

(3.11) 
(3.12) 



where the second term in the brackets of ( |3.11| ) contributes to the subleading terms denoted 
by ellipses in (|3.12|) . 

Secondly, considering case ii), we find that differentiation gives two terms which will 
cancel in the leading order, 



lim Jim A (a)< (P',P) 



ig 



57P 



Jo dU J-i *^ 74 ~ V ^ ' $ 



Sign(0-^V(lX) + 26(OV(lX) 



(3.13) 
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(a) (b) (c) 

FIG. 2. The physical radiative corrections to the quark-gluon vertex; a) A^ , the abelian 
vertex, b) A^ b \ the tri-gluon vertex and c) A)} , the triangular vertex. 

The first term is identical to that evaluated for case i) above. With the same approximation, 
the second term is, 

" t dl n J 2 a , , - ^ (3.14) 

4tt 2 7o Z 3 + fmf>| 12tt 2 7rmf>| V 7 

The two leading contributions cancel and in this ordering of limits the vertex is finite. 

Now we can see how the paradox is resolved. The spatial abelian vertex considered with 
the ordering of the limits in case ii) is finite, in agreement with the second part of the 
identity (O). 



B. The BRST Identity. 

The BRST identity is a generalisation of the Ward-Takahashi identities for non-abelian 
gauge theory obtained through the BRST transformations. The BRST version of the Taka- 
hashi identity can be written as, 

(P' _ Aj^P', P) = gT l f {Z{P') - S(P)) + R l (P', P). (3.15) 

The physical quark-gluon vertices = Aj^+A'w+Ajyy are represented in Fig. |2|. The non- 
physical ghost-quark vertices induced by the BRST transformation, R l = Pj^ + Pjw +R l ( c \ 
are represented in Fig. |[ They vanish for on-shell Minkowski momenta P and P' at fi = 0. 

The nontrivial part of the BRST identity ( |3.15D is in the dressing of the gluon lines of 
Figs. [1], ^] and |3] by HDL and the inclusion of Fig. |2]c. The order of the perturbative expansion 
is mixed up without offsetting the simple form of the identity. The detailed derivation of 
( |3.15| ) is given in Appendix A. 

Setting p ' = p on the Fermi level, the BRST identity (|3.15| ) implies that 



lim 



V — V 



=p, p=fj. 

,3 /| ;3 



It follows from the discussions of the previous subsection that A{(P',P) = A7 j (P',P) + 
Af\P',P)+A l } c \p',P) with 
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(a) (b) (c) 

FIG. 3. The non-physical ghost diagrams generated by the BRST transformations; representing 
o) R l{a \ b) R l( & and c) R 1 ^. The open circles denote non-physical vertices generated by BRST. 



kf\P\ P)=gT\ 



a 



ad 



+ c t 



ig 2 , 41 
7^74 In — 



Txm D \v 



(3.17) 



It remains to find the logarithmic terms from Figs. 0b, 0c or |3| to reconcile the BRST identity 
Q3.15Q to the leading order of the infrared logarithms. 
For p ' = p and p = //, we have 

AjfV, P) = ig"f alh T^T b f - £ J 7^(2^ -v'- v)V% u n )^ - Ui)liS{P + Lfri 



^gC ad T l f [A ( r(P', P) + Ai 6)> (P', P)} 



(3.18) 



where the superscripts specify the contributions from loop momentum inside and outside 
the infrared region, |a>| < uj c and I < l c . The inside contribution can be approximated by 



where 



and 



Af }< (P',P) 



47T 3 



dll 



1 



-I i(uj + v) — £ 



2tt 3 



[ihlA + hl-p), 



I x = dll dcu toV 2 (I, u) tan" 1 

JO J —U) r 



I 



LO + V 



(3.19) 



(3.20) 



dll / dLUuV 2 (l,u) 



0J + V 

I 



tan 



uj + v 



(3.21) 



In the limit v — > 0, the integrand of both integrals I\ and I2 are positive and can be bounded 
by letting u c — > 00 and changing the integration variables from to (l,x = u/l). We 

have 



r°° XI 
\h\ < / dx— 
Jo 1 



x x tan 1 - 



+ x 



2\2 



hi 



l 2 c (l + x 2 ) 
m 2 D f M (x) 



(3.22) 



and 
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\h\ < / dx 



oo x (x — tan 1 - 



o (1 + x 



2\2 



m 2 D f M {x) 



(3.23) 



Both integrals are convergent and hence does not contribute to the infrared logarithm. 
It is also straightforward to verify that the BRST generated diagrams in Fig. |3] do not 
display any logarithmic behavior in the limit v — > 0. Thus the only candidate left over is 
the diagram in Fig. 0c corresponding to gluon insertion on a HDL. 

Though formidable as it looks, evaluation of Fig. ^|c can be simplified with the aid of a 
Ward type identity which relates the derivative of the gluon self-energy and the tri-gluon 
vertex with three external gluon lines. Again the answer is sensitive to the relative order of 
the limits v' — > v and p' — > p. In appendix B, we shall demonstrate that 

lim \imA l } c \P',P) = ^7-5^74 In -^2— + . . . , (3.24) 
v'->vp'-^p 2Aty z Tim D \i'\ 

which completes the BRST identity (|3.16|) to the leading logarithm level. The other order 
of the limits, linv-^ \m\pi P^j (P', P) is infrared finite. 

To summarise, we have shown how the BRST identity works in a quark-gluon plasma at 
high density in terms of Feynman diagrams. The derivation of the identity is quite general 
and that it works at high density is not surprising, but with this presentation we hope that 
the mystery that shrouds this topic may be lifted. With the incorporation of HDL, the 
BRST identity is no longer satisfied order by order. The payment for using HDL is that 
orders of perturbation theory become mixed up. 



IV. COLOR SUPERCONDUCTIVITY 

Perturbative QCD has been applied successfully toward the study of color supercon- 
ductivity at high baryon densities. In this regime, single gluon exchange dominates the 
pairing interaction, and screening plays an important role in the non-BCS behavior of color 
superconductivity [|l^,|l^-[T7| . In Refs. [fI8| , pO|l , the superconducting pairing temperature of 



a dense quark-gluon plasma was investigated by means of a Dyson-Schwinger approach to 
the pairing interaction. The resulting problem was reduced to one of finding the smallest 
eigenvalue, A, of the Fredholm equation 

A 2 f°° 

fs' lS ' 2 {n'\p) = — / d P K s' 1 s' 2 ,s 1 s 2 {n',n\p',p)f SlS2 {n\p), (4.1) 
with the condition A 2 = 1 yielding the critical temperature. The kernel is given by 

K s ' iS ' 2iSiS2 (n',n\p',p) = — 22 ls' lS > 2 , s '{ s >>(n',n\p',p)S s >> Sl (n\p)S s >> S2 (-n\p), (4.2) 

and consists of the s-wave components of the two particle irreducible amplitude for the 
scattering of two quarks in their color antisymmetric channel with zero total energy and 
momentum, 'j s ' 1 s' 2< s 1 s 2 (n' , n\p',p), and the full quark propagator S s > s (n\p). The initial energies 
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T(n',n\p',p) = 




S(n\p) = 



FIG. 4. The diagrammatic expansion of the two particle irreducible vertex T s i s / s s (n',n\p',p) 
to order g A and the quark self-energy S s i s (n\p) to order g 2 . 



of the two quarks are ±iz/ n and the final ones are ±izv with v n = (n + |)e. The initial 
momenta of the two quarks are ±p and the final ones are ±p'. The diagrammatic expansion 
°f 7s' 1 s' 2 ,s 1 s 2 ( n ' , n\p',p) to order g 4 and S s i a (n\p) to order g 2 is shown in Fig. |4[ 

Collecting previous results, the perturbative expansion of the least eigenvalue reads 

m 

,2 



1 

A* 



9 



24tt 5 



— log 2 - + — ( 7 + log 2) log - + 0(1) 



7T Z 



7T Z 



24vr 2 



1 + 



' 4(7r 2 + 4). 
C f _ A ' log 



7T J 



3 2 + o(iog 2 -; 



(4.3) 



where the leading 0(g 2 ) term stems from the first diagram of Fig. |] with a bare quark 
propagator. Relative to this leading term, the radiative corrections and the two gluon 
exchange appear to be suppressed by g 2 as is the case with the remaining diagrams of Fig. £|, 
but may not be so because of the infra-red logarithm, each counted as g~ l for T ~ Tc in 
accordance with ( |1.2j ). The radiative correction to the quark propagator is such an example 
The logarithm of the self-energy, contained in the second line of 
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• 3|), gives rise to a 

significant contribution to the prefactor JL8|. Though the radiative correction to the vertex 
function is liable to such a logarithm according to the BRST identity, this does not happen 
in the energy momentum region \p — p\ ~ \v\, \p' — p\ ~ \u'\ and \p — p'\ ~ [k\v — f'Ds, 
where the main contribution to the kernel (|4.2|) comes from; this is indicated by the absence 
of the logarithm in the limit v' — > v followed by p' — > p of the vertex function. In what 
follows, we shall demonstrate this point via an explicit evaluation of the contribution of the 
abelian vertex function to the partial wave amplitude. 



Consider the abelian vertex Fig. |^a, with p 
contribution is given by 



p 



p, q = 2p sin | ~ p9. The infra-red 



A^(P',P) =u(P')kf\P\P)u{P) 



9 2 Pj 



9 2 Pj 



d 3 t 
(2^)3 

(2^)3 



duj 

—u) c 27T 

duj 
c 2% 



p-p' - {p-i){p' ■/) 

p ■ p' — (p ■ I) {p' ■ I) 



(4.4) 



iAv + ( 



1 



where Au = i/ — u, £ = \p+l \ —p and ^' = \p'+l \ —p with |^| < / and |^'| < I. It follows from 
the discussions of the previous sections that the sensitive region of the integration variables 
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which is responsible to the non-Fermi liquid logarithm corresponds to the singularities of the 
fractions inside the bracket. Therefore one of £ and £' must be kept small in the sensitive 
region. If there were an infra-red logarithm, it would come from 



A (a) (p' p)=i — r 

" 1 ' } Ji<i c (2^J- w 



d 3 l pc du T>(l,u) 



(2vr) 3 J- Wc 2vr iAu - f ' + £ 



9(rjl-\^\) 9( V l-W\) 



(4.5) 



with 77 1. Transforming the integration variables from I to I, £ and we have 



d 3 l = ?-ldl (4.6) 



where the Jacobian is J = \l • p x p'\ ~ fi 2 ^Jl 2 9 2 — (£ — £') 2 with the approximation that £ 
or £' <C /. Introducing 

P(Z,w) = r dw'V(l,u'), (4.7) 



and carrying out the integration over £ and we obtain 

Note that if p ' — ► p first, we have a complete exposure of Az/ in the denominator, 

(4-9) 

The integration will give rise to z/'ln|z/'| — i/ln|i/|, which in the limit v' — > v produces 
the infra-red logarithm. But here, with 9 ~ |z/ — i/|s and / ~ |z/|3, such a singularity is 
suppressed through the l 2 9 2 term inside of the square root. Indeed, if we insert Atf>(P',P) 
into the partial wave integration, we find the corresponding contribution 

™<"» = 5? jf *»W " * £* ?pirnz (410) 

7/ii( z// ) ^) ma y be bounded by dropping Az/ inside the square root. Then the integration 
over 9 decouples from that over / and u, i.e. 

\1ir(v',v)\<£- 2 IJ, (4.11) 

where 

/= f C d9V(fi9,u' - v) (4.12) 

JO 

and 
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J = [° dl I V duV(l,u). (4.13) 

JO J-u' 

It follows from the properties of the function D(l,u) that I and J are bounded from above 
by 

/< -^(«|i/-i/|)-s (4.14) 



and 



J < -?L K -*\\j/\* - Sign(z/V)|z/|i|, (4.15) 
v3 

where k = \m 2 D . Combining / and J, we see that 7t\r( z/ ' ) v) is nonsingular in the limit 
v' —>■ and v — > along any path in the (z/, ^)-plane. 

It is important to note that this result for 7/i?(z/', v) only pertains to any possible ad- 
ditional infra-red enhancement arising from the radiative quark-gluon vertex A l J> a >. For the 
complete partial wave amplitude, 7(1/, u), the collinear magnetic gluon exchange logarithm, 
already present at tree level, i.e. 

7tree (z/,z/)^ln^, (4.16) 

maintains its presence at the radiative level. Indeed, based on numerical evaluation of the 
partial wave amplitude, 7abeiian(z // , v), we have confirmed that only this expected collinear 
logarithm is present. We have also evaluated the infra-red contribution 7/_r(z/', v) numeri- 
cally, and the result supports the above analytic arguments. 

Though our conclusion that the vertex function does not contribute to the pre- 
exponential factor agrees with that made in |T5|, the arguments used in |TS[ to justify this 



conclusion merit further consideration. In particular, the formula for the vertex function, 



taken from Ref . |y] , is not applicable for the infra-red contribution at a large chemical poten- 
tial in comparison with the temperature. This can be judged by the absence of the infra-red 
logarithm from their vertex function in any order of the limit of zero energy-momentum 
transfer; this absence contradicts the BRST identity as discussed above. In fact, only the 
expressions for diagrams with internal fermion lines only can be carried over from the high 
temperature region to the large chemical potential region, as is the case with the gluon self- 
energy functions ( |2.5| ) and (|2.6|) . For diagrams with internal gluon lines, the infra-red region 
makes significant contributions, leading to effects such as the non-Fermi liquid behaviour of 
the quark self-energy and vertex functions, which has been completely ignored by the Hard 



thermal loop approximation employed in |£1 



By careful examination of the radiative corrections to the quark self-energy and vertex 
functions, we have reconciled the non-Fermi liquid behavior in the dense plasma with the 
BRST identity. The incorporation of HDL, and the resulting resummation in the gluon 
propagator, leads to a mixing of orders in the perturbative expansion. Hence proof of 
BRST involves combining diagrams of different loop order, as seen in Fig. |^. An important 
consequence of this result for color superconductivity is the verification that there are no 
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additional infra-red logarithms accompanying the radiative correction to the vertex function. 
This strengthens our previous result that the only radiative correction to the determination 
of To comes from the quark self-energy, and suggests that the pre-exponential factor of ( |1.3| ) 
is in fact exact to leading order in g. 
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APPENDIX A: 

In this appendix, we shall prove the BRST identity, (|3.15p , relating the self-energy, vertex 
and ghost diagrams of Figs. [TJ, ^| and ^ in the presence of hard dense loops. 
Using the standard trick, 

(P 7 - PyS{P' + L) lfl S(P + L)= S(P' + L) - S(P + L), (Al) 

we may trivially relate the abelian vertex Fig. ^a with self-energy Fig. |1]. 

(P ' _ p^Af^P', P) = gT) (\ - Q±) (E(P') - £(P)) . (A2) 

which, apart from the group theoretic factors, is nothing but the Takahashi identity of QED 
and is independent of the form of the gluon propagator. However, for non-abelian gauge 
theories, the group coefficients do not match; cancellation of the extra term must result from 
the additional vertices. 

In QCD there is a second physical process at 0(g 3 ) in perturbation theory, formed with 
the tri-gluon vertex —if lmn T fl xp, as shown in Fig. 0b. We now turn to this diagram to see 
how it may cancel the extra terms. It is straightforward to write down an expression, 

Af\p>,P) = f™T™TfjY: J ^V vX {L){-t^{L,L-Q)V^{L) lv S(P + L) lv ,, 

(A3) 

where if lmn jj l j 1 ^ = _<Z|d j'j. However, this expression may be simplified when contracted 
with (P' — P)p with the aid of the identity 

(p> - p)^ A (p')rWP', p)v pv {p) = i {k2(p', p) + vW(p>, p) + v®(P', p)} , (A4) 

where 
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K2(P', P) = i [DMP) - VMP')] , (A5) 

v^(p', p) = v v , x {p') [u Xp (p') - n v (P)] v pv {P\ 

V®{P', P) = A{P')P^V Xv {P) - V U , X {P')P X P U A{P). 

n^P) is the HDL diagram which satisfies P p Ii^ u {P) = and A(P) = —i/p 2 is the ghost 
propagator. Since II is itself of 0(g 2 ), we see here that the price one pays for incorporating 
HDL in the gluon propagator is the mixing of orders in perturbation theory. To prove ( |A4|) , 
we start with the bare gluon propagator, 

D,AP) = p 1 2 [^ + (a-l)^fl (A6) 

and the identity 

(P 1 - P)»(-i)T wX (P', P) = (P 2 - P' 2 )5 pX + P' p P' x - P P P X . (A7) 

Sandwiching ( |A7|) between D(P') and D(P), we find 

- i{P' - P) p D a , p {P')T, pX {P\ P)D a {P) (A8) 
= -i[D a , a {P') - D a , a {P)) + A(P')P^P' p D pa (P) - A(P)D ap (P')P p P a . 

The HDL-dressed gluon propagator is related to the bare propagator via the Dyson- 
Schwinger equation, 

V pu (P) = D pv (P) - iD pp (P)Tl pX (P)V Xu (P) = D pu (P) - tV pp (P)U pX (P)D Xl/ (P). (A9) 



It follows that -i(P' - P) p V v , x (P')T pXp (P' ', P)V pu (P) can be obtained by sandwiching (|A8 



between 8 v i a i — %D v i p (P')Upi a i (P 1 ) on the left and 8 au — iH al 3(P)Vf 3u (P) on the right. The 
expression is then simplified by the 4-dimensional transversality of the self-energy matrix 
n(p), and we end up with (|A4|) and (|A5|). 

Now we look at the contribution due to V^ l \ 

'^f T fJ E / -Q)- VML)} 1,S(P + L) lv , = ig^-T) (E(P') - E(P)) . 

(A10) 

This expression, the origins of which are purely non-abelian in nature, exactly cancels the 
extra terms induced in ( |A2| ) . However, the tri-gluon vertex also induces a number of extra 
terms which we shall now consider in turn. 

The appearance of the ghost propagators in suggests that these extra terms will be 
cancelled by the non-physical ghost-quark vertices generated by the BRST transformations. 
Indeed, we find that, when grown from a fermion line, the ghost terms contribute, 

i?™TJTf^ £ J j0r 3 V^(L, L - Q) lv ,S{L + P) 7 „, (All) 

which exactly cancels the diagrams of Fig. |3|. 

17 



FIG. 5. Quark loop with three external gluons, rjjjjjj. 



The remaining term, V^ 2 \ is of 0(g 5 ), two orders higher in perturbation theory, and 
contributes, 

i?™T?Tf^ £ J 703^a(£ - Q) [K Xp {L) - Tl\ p (L - Q)\ V pv (L)^S(L + P) lu . (A12) 

In the absence of HDL this term does not appear and the BRST identity is satisfied order by 
order in perturbation theory. Although with the inclusion of HDL the ordering has become 
mixed up, the identity must remain. To see how this contribution is cancelled we study the 
triangular vertex shown in Fig. |2]c. 

We shall first look at one of the two loops that form the triangular vertex, namely the 
quark loop with three external gluons, as shown in Fig. With three identical vertices 
there are two possible orderings for this diagram. Considering both orderings, we write this 
vertex correction as 

f X( p/ > P ) = ^E / ^Tr [t\T?T? 1x S{L + P - P') ltl S{L) lp S{L + P) (A13) 

+ T l f TpT? lx S(-L - P) lp S(-L) llM S(-L -P + P') 
Contracting one leg of the triangular vertex with (P 1 — P) M we find, 

(p' - Pff X(p', p) = igf lmn [IMP') - imp)] , (Ai4) 

where, as discussed, n is the vacuum polarization diagram. Therefore, connecting the two 
free legs to a fermion line, we find that 



(p' -pykf\p',p) = 

-fmn T p T n^ £ J ^-V^L - Q) [U Xp (L - Q) - U Xp (L)} V pv {L) lv ,S{L + P) 7 „, 

(A15) 

which cancels the remainder in ( |A12 ). 

At this point, we note that at nonzero chemical potential the triangular diagram, ( |A13| ), 
does not contain exclusively the term proportional to f lmn j and is nonvanishing even in 
QED because of the breakdown of Furry's theorem by the Fermi sea. On the other hand, 
the identity ( gig) remains valid rigorously and there is no contribution from the triangular 
diagram to the Takahashi identity of QED. Furthermore, for low excitations near the Fermi 
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level, the approximate particle-hole symmetry renders the triangular diagram dominated by 
the term proportional to f lmn . 

Before concluding this appendix, we shall relate the particular BRST identity (|3.15 ) 
to the master BRST identity as given in Ref. [25|. Let T(A, x, X, c i c) be the generating 



functional of proper vertex functions with A, Xi x, c an d c the quantum mechanical average 



of the gauge potential, Vu, quark fields ip, ip and the ghost fields 



i.e. 



reads 



(ip(x)), x( x ) — W^))) c ( x ) — (0( :r )) an d c ( x ) = The master BRST identity 



d x 



5T 



SA^x 



<8VJx)) + 



where 



(A16) 



SV. 



dx, 



+ gf 



■in \/'" () '' 



5^ = -iT l (f> l ip, 
Sip = -iT l <\> l i\), 

= _^mn^m^n^ 



(A17) 

(A18) 
(A19) 

(A20) 



are the BRST variations of the field components. The expansion of the term xx c i n ( A16 ) 
to the order g 3 and with the bare gluon propagators replaced by the dressed ones afterwards 
yield the identity ( |3.15|) . Unlike an abelian gauge theory, the ghosts couple to other fields of 
the theory. The expectation of the nonlinear term of the BRST variations gives rise to the 
additional terms R a (p',p) with R 1 ^ from the second term of ( |A17|) , R 1 ^ from QA18| ) and 
R 1 ^ from (TAT9D . 



APPENDIX B: 

In this appendix, we shall evaluate the infra-red contribution of the diagram in Fig. ||]c, 
which we denote by K l !f\P' , P)ir with P = (p, v) and P' = (p + q,v + Av). Then Q = 
P' — P = (q,Au), and both q and Av are soft. The calculation is greatly simplified with 
the aid of the identity ( |A14f) for // = 4 in the limit q — > followed by Av — > and for \x = j 
in the limit Av — > followed by q — > 0. 

(i) The triangular vertex in the limit 



We start with 
kf\P\P) m = -g" 



d 3 l 



l<l c (27r) 3 J-ui c 2tt 

x V(\l- q\,u- Av)V{l,u) 



lim limAi (c) (P',PW 

Ai/^O <f->0 



^T a T b [-iT l C n/ (L,L-Q)]^ 



(Bl) 



P + l 



In 



{I ~ q) m >{l - q), 
\l- q\ 2 



{$n'n — ln.dn)i (B2) 
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where T>(l,u) is given by (|2.8|) and rj^" by Fig. §. Note that we have used the continuum 
approximation for the Matsubara sum. It follows from ( |A14j ) that 

lim hmfg(L,L - Q) = igf^Ik^L). (B3) 

Av— >0 q^O J OOJ 



Therefore 



lim limA? \P', P) IR = \gC ad T l A{P), (B4) 



with 



v ' Ji<i c (2vr) 3 J-u c duo v ' i(u + v) - i K ' 

ig 2 /<■.„ r , da M i ft .e^ 1 



8tt 3 7 J-^ duo [I 2 + UJ 2 + a M (l,u)] 2 J-i ^' I 2 1 *'i(w + v)-£' 
where p = \i and £ = \p + I \ — /i. Carrying out the £ integration, we find 
Arp s ig 2 & dl r da M Fjy^u) 

A(p) = 2Wo 7U^-fo[r+u?+*"{i,u,)r (B6) 

with F(v, p,;l,u) given by (P-14 ) The discontinuity of the inverse tangent corresponds to 
uj ~ — v and the /-integration is dominated at I ~ (kuj)3 ~ [kv)3. We end up with 



g 



Att 3 Jo I J-^ du [I 2 + a M (l,u)} 2 

2 Ac 

74 / dlT>(l, —v) + terms regular as v — > 
Jo 



4tt 2 

^ 2 , i 



74 log — ^— + terms regular as v — > 0. (B7) 

127T 2 K\U\ 

(ii) The triangular vertex in the limit 

lim lim A l } c) (P',P) IR . (B8) 

q— >0 Av— >0 

Here A l } c) (P',P) is given by ( [B"2] ) with the replacement T l £ff n , — > Y l ™J n ,. It then follows from 
the identity (gig ) that 

limlimf;t(L,L- Q) = igf lab ^-U mn (L). (B9) 



Therefore 



hm Hm, Af \P', P) IR = ^gC ad T%A>(P), (BIO) 
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with 
A'(P) = 



fr 



8tt 3 

-.2 



duj 



da 



9 



L dl L^~ di \i 2 + uj 2 + 

df M i 



i r l t 2 l 

—^J- / dt^-i — l-p)- ■ r 

■ a M (I, u)\ 2 J-i V i{uj + v) — 



47T 3 



d fdl f 

JO J —LJ C 



dcu- 



dl 



l 2 + u 2 + m 2 D f M ^ 



x j 74 
+i 



i 

-1 H : — tan 1 



U) + v ) 

\1\ I, , i ,A3 



7-p . (Bll) 



The discontinuity of the inverse tangent function at uj + v = is now smeared by the factor 
u) + v. This integral converges in the limit v — > 0. 
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